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Abstract 

We examine a family of 3-point mappings that include mappings solvable through lin- 
earization. The different origins of mappings of this type are examined: projective 
equations and Gambier systems. The integrable cases are obtained through the appli- 
cation of the singularity confinement criterion and are explicitly integrated. 



1. Introduction 



Integrability is a term far too general. Although the idea is simple and related 
to the integration of a differential system, integrability assumes various forms. Here 
we shall attempt neither a rigorous definition nor an exhaustive description of all the 
disguises of integrability. In order to fix the ideas we shall just present three most 
common types of integrability, which suffice in order to explain the properties of the 
majority of integrable systems [1]. These three types are: 

- Reduction to quadrature through the existence of the adequate number of inte- 
grals of motion. 

- Reduction to linear differential systems through a set of local transformations. 

- Integration through 1ST techniques. This last case is mediated by the existence 
of a Lax pair (a linear system the compatibility of which is the nonlinear equation 
under integration) which allows the reduction of the nonlinear equation to a linear 
integrodifferential one. The above notions can be extended mutadis mutandis to the 
domain of discrete systems. 

This paper will focus on the second type of integrability, usually refered to as 
linearizability. The prototype of the linearizable equations is the Riccati. In differential 
form this equation writes: 

w' = aw 2 + (3w + 7 (1-1) 
and the transformation w = P/Q (Cole-Hopf) reduces it to the linear system: 

P' = f3P + 1 Q (1.2) 

Q' = -aP 

Similarly, the discrete Riccati equation, which assumes the form of a homographic map- 
ping: 

_ ax + /3 

x = — (1.3) 

where x stands for x n , x for x n+ i (and, of course, x for x n _i), can be linearized through 
a Cole-Hopf transformation. Putting x = P/Q we obtain readily: 

P = aP + (3Q (1.4) 

Q = 7 P + 5Q 

The extension of the Riccati to higher orders can be and has been obtained [2]. The 
simplest linearizable system at N dimensions is the projective Riccati which assumes 
the form: 

w 'u = q m + ^ b^Wy + W/j c u w v with fj, = 1 , . . . , N (1.5) 

V V 

In two dimensions the projective Riccati system can be cast into a second order equation 
without loss of generality: 

w " = -3 WW ' - w z + q ( t )(w' + w 2 ) (1.6) 



The discrete analog of the projective Riccati does exist and is studied in detail in [3]. 
The corresponding form is: 




(1.7) 



Again in two dimensions (and only in this case) the discrete projective Riccati can be 
written as a second-order mapping for a variable without any simplifying assumptions. 
One is thus led to a mapping of the form 



which was first introduced in [4]. The coefficients a, ft, ... ,9 are not totally free. Al- 
though the linearizability constraints have been obtained in [4] , the study of mappings of 
the form (1.8) was not complete. In the present work we intend to present its exhaustive 
study from the point of view of integrability in general. 

Another point must be brought to attention here. In the continuous case the study 
of second order equations has revealed the relation of the linearizable equations (1.6) to 
the Gambier equation [5]. The latter is obtained as a system of two coupled Riccati in 
cascade 



where n is an integer. The linearizable equation (1.6) is obtained from (1.9) for n = 1 
and a = — 1, c = and a = 0. The discrete analog of the Gambier mapping was 
introduced in [6] and in full generality in [7]. However the relation of the linearizable 
mapping to the discrete Gambier system has not been studied in full detail. In what 
follows we shall fill this gap by presenting the reduction of the Gambier mapping to the 
linearizable one. 

In the next section we shall analyse (1.8) and all its reduced forms and isolate the 
integrable ones through the use of the singularity confinement criterion. The integration 
of all integrable cases will be given in detail and in particular the linearization through 
a projective system. The third section will be devoted to the study of the reduction 
of a Gambier-type (coupled Riccati) system to a linearizable one together with the 
investigation of the integrable cases and their integration. 

2. Linearizable mappings as projective systems 

In [4] we have introduced a projective system as a way to linearize a second-order 
mapping. (The general theory of discrete projective systems has been recently presented 
in [3]). In this older work of ours we have focused on a three-point mapping that can 
be obtained from a 3 x 3 projective system. As a matter of fact, this is the only case 



axxx + [3xx + 'yxx + 5xx + ex + (x + rjx + 9 = 



(1.8) 



y = -y 2 + qy + c 



(1.9) 



w' = aw 2 + nyw + a 



where the projective equation can be converted to a single, one component, mapping 
without any simplifying assumptions. The main idea was to consider the system: 
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(2.1) 



and conversely 
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m 13 
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m 33 



(2.2) 



where the matrix M is obviously related to the matrix P through M = P _1 . Intro- 
ducing the variable x = u/v and the auxiliary y = w/v we can rewrite (2.1) and (2.2) 
as 

_ Pnx + p 12 + p l3 y 

x = (2.3) 

p 21 x +p 22 +p 23 y 

_ miio; + 77112 + m 13 y 
~ m 2 \x 4- m 2 2 4- m 23 y 

(Since the 7773^, p 3 i do not appear in (2.3), (2.4) and we can simplify M,P by taking 
77733 = P33 = 1 and 77731 = 7731 = 77732 = 7732 = 0). Finally eliminating y between (2.3) 
and (2.4) we obtain the mapping: 



axxx + (3xx + "fxx + dxx + ex + (x + rjx + 9 = 



(2.5) 



where the a, (3, ... 9 are related to the m, p's. 

Equation (2.5) will be the starting point of the present study. Our question will 
be when is an equation of this form integrable? (Clearly the relation to the projective 
system works only for a particular choice of the parameters). In order to investigate 
the integrability of (2.5) we shall use the singularity confinement approach that was 
introduced in [8]. The first question is what are the singularities of (2.5)? Given the 
form of (2.5) it is clear that diverging x does not lead to any difficulty. However, another 
(subtler) difficulty arises whenever x n _|_i is defined independently of x n -\. In this case 
the mapping "loses one degree of freedom" . Thus the singularity condition is 

dxn+i = Q 
dx n -i 

which leads to : 

{ax 4- 6) (ex + 9) = (fix 4- C){l x + v) (2-6) 

Equation (2.6) is the condition for the appearence of a singularity. Given the invariance 
of (2.5) under homographic transformations it is clear that one can use them in order 
to simplify (2.6). Several choices exist but the one we shall make here is to choose the 
roots of (2.6) so as to be equal to and 00, unless of course (2.6) has two equal roots, 



in which case we bring them both to 0. Let us examine first the distinct root case. For 
the roots of (2.6) to be and oo we must have: 

ae = 07 (2.7) 

66 = & 

The generic mapping of the form (2.5) has a6 ^ and we can take a = 6 = 1 (by the 
appropriate scaling of x and a division) . We have thus, 

xxx + (3xx + ■yxx + (rjxx + (3^x + (x + r\x + 1 = (2.8) 

Nongeneric cases do exist as well and we shall examine them in detail. Thus let us 
first assume that a = in which case (3^ = and we can decide that (3 = without loss 
of generality {(3 and 7 are interchanged if one reverses the direction of the evolution). 
The mapping then becomes: (7a; + e)x + ((x + l)(r]x + 1) = and we must have ( 7^ 0, 
lest the mapping become a two-point one, in which case a scaling of x can bring its 
value to 1. We have thus the general form: 

xi^x + e) + (x + l)(r]x + 1) = (2.9) 

Clearly we cannot take both r\ = 0, 7 = 0, neither e = 0, 7 = 0. So, in order to obtain 
the reduced forms of (2.9) we assume first 1] = leading to 

x + 1 + x(~/x + e) = (2.10) 

and then rj 7^ with two possible combinations either 7 = 0: 

(x + l)(r]x + 1) +ex = (2.11) 

or e = 

(x+ 1)(t]x + 1) +-fxx = (2.12) 

Going back to mapping (2.5) we can obtain a further reduced form by assuming a = 
6 = 0, in which case (3^ = and (?] = 0. As previously, we can assume that ,9 = but 
this fixes the direction of evolution so the two choices ( = and 1] = are distinct. We 
have thus two reduced mappings: 

r yxx + 5xx + ex + (x = (2-13) 

and 

jxx + Sxx + ex + r]x = (2-14) 

A further reduction can be obtained if we assume that (3 = 7 = 0, in which case we can 
freely choose ( = and find: 

Sxx + ex + t]x = (2-15) 



Mappings (2.8) to (2.15) are all that can be obtained from (2.5) with the assumptions 
of the distincts roots for (2.6) up to homographic transformations. 

Next we turn to the case where (2.6) possesses a double root equal to 0. The 
constraints in this case read: 

89 = Cv (2-16) 
a + Se = (hi + Cl 
The generic case corresponds to 9 = e = 1, leading to : 

axxx + j3xx + "fxx_ + (rjxx + x + Qx + r]x + 1 = (2-17) 

with a = f3r] + (7 — 5. Nongeneric cases exist also and we start by considering 9 = 0. 
This leads to £77 = and we choose £ = 0. One further constraint must be satisfied in 
this case: 5e = (5r\. Assuming e = 1, we have 8 = (5r\ and obtain: 

axxx + (3xx + "fxx + ftrfxx + x + rjx = (2-18) 

If we take e — then two choices exist: rj = or f3 = 0. The first leads to: 

axxx + (3xx + "fxx + 8xx = (2-19) 

while the second gives: 

axxx + "fxx + 8xx + rjx = (2.20) 

We can immediately point out that the mappings (2.19) and (2.20) are trivial. The first 
reduces to a affine one under the transformations x — > 1/x, while in the second the term 
a; is a common factor and the mapping is in fact a two-point one. The mappings (2.17) 
to (2.20) complete the reductions of (2.5) in the case of a double root of (2.6). 

In order to investigate the integrability of the mappings obtained above we shall 
apply the singularity confinement criterion as previously explained. We start with the 
generic mapping (2.8). Here the singularities are by construction and 00. Following 
the results of [4] we require confinement in just one step. These leads to the condition 
(3 = ( = 0. We obtain thus the mapping: 

xxx + 'jxx + r]x + l = (2-21) 

or, solving for x : 

W= _ 1+ !K±1 (2.22) 

where 7 and r\ are free. We expect (2.22) to be integrable. This is indeed the case: we 
can show that (2.22) can be obtained from the projective system (2.1), (2.2) provided 
we take 

-(79 + 1) QQ 1 
P= I q I (2.23) 

1 



and M = P _1 provided q is some solution of the equation qqq + qqr\ + 57 + 1 = 0. 
Then we can use the projective system (2.1) and (2.2) in order to construct the general 
solutions of (2.21). (Let us point out here that the equation for q is really different from 
(2.21) due to the presence of 7 rather than 7). 

We turn now to the remaining mappings starting with (2.9). For the application 
of the singularity criterion we must find the values of x for which the mapping loses one 
degree of freedom. For mapping (2.9) this happens (by construction) when x is or 
00, but another source of singularity exists, when 72; + e = 0. The investigation of the 
singularity x = leads to the confinement condition e = rj = 1. Proceeding now to the 
examination of the singularity starting from 'jx + 1 = we find that there is no way to 
confine it: the sequence is {— 1/7, x, —1, /, 00, 00, . . .} where / is some finite value. We 
conclude that (2.9) is not integrable. 

Next we examine mapping (2.10): 

x = -1 - x(~/x + e) (2.24) 

If x=0 we find x = — 1 and all the subsequent values are independent of the initial data: 
there is no way for the mapping (2.24) to recover the lost degree of freedom and thus 
we conclude that this mapping is not integrable. Next we turn to (2.11): 

x = -l ^— (2.25) 

TjX + 1 J 

and again start with x = 0. We obtain x = —1, x = — 1 + e and then x may become 
indeterminate, in the form 0/0, provided e = 77 = 1. A careful analysis of this case 
shows that the singularity is periodic. Translating x we find a mapping of the form: 

xx + x-l = (2.26) 

which must be integrable. It turns out that (2.26) is indeed integrable. Its integral can 
be readily obtained: 

K= z + ^ + J-- x -x- 2 (- + -) (2.27) 



rp ry> 'T* ' T* \ nf If 

and (2.27) can be integrated in terms of elliptic functions. Mapping (2.26) is a member 
of the QRT family of integrable mappings [9] and can be easily shown to be the only 
member of this family included in the parametrization (2.5). Mapping (2.12): 

- = _ x _jx^_ 

T]X + 1 J 

has two kinds of singularities induced either by x = or r\x_ + 1 = 0. The latter leads 
to an unconfined sequence x = 00, x = 00, etc and thus we conclude that (2.28) is not 
integrable. We turn now to mapping (2.13): 

_ = _ x jx^+e 
Sx + C 



A singularity starting with 5x+( = leads to an unconfined sequence x = oo,x = oo,. . . 
unless 7 = 0. However even if 7 = we obtain a sequence 00, 00, finite, 0, 0, . . . again 
without recovering the lost degree of freedom. Thus this mapping too is not integrable. 
Next we consider mapping (2.14). If 5 7^ we can put 5 = 1 by division (if (5 = the 
mapping becomes a two point one): 

ax 

x = -(- fX + r)) (2.30) 

The singularity of (2.30) occurs when x = and we obtain succesively x = — 1 and 
x = 00. The condition for the subsequent x's to be finite is just e = 777 or e = 0. The 
latter is trivial because the mapping becomes a two-point one and we consider only the 
first condition. We take rj 7^ 0, in which case a scaling can bring its value to rj = 1, 
resulting to e = 7. Thus, based on singularity confinement the mapping: 

x(<yx + x + 1) +7X = (2.31) 

should be integable. This is indeed the case. Its integration can be obtained in a 
straightforward way. We find as the integral of (2.31) the quantity: 

(x + l)(x + l) 
ax 

where a is related to 7 through 7 = —a/ a. Thus (2.31) is the discrete derivative of 
a homographic mapping. Once (2.32) is obtained the complete integration of (2.31) is 
elementary. 

Finally, to make the long story short, we can just give without detailed proof 
the result that none of the mappings (2.15), (2.17) and (2.18) satisfies the singularity 
confinement criterion. Thus we expect all of them to be nonintegrable. 

Before closing this section let us present the continuous limits of the two lineariz- 
able equations we have identified above. For (2.21) we put x = —1 + vw, 7 = 3 + v 2 p, 
rj = 7 + v 3 q and we obtain at the limit v — > the equation: 

w " = 3ww' - w 3 +pw + q (2.33) 

This is equation ^6 in the Painleve/Gambier classification [10] (in noncanonical form) 
and precisely the one that can be obtained from a N = 2 projective Riccati system. For 
equation (2.31) we put x = 1 + vw, 7 = — 1 — u 3 q and obtain, at the limit v — > 0: 



(2.32) 



w " = ww > + 2q (2.34) 

This is equation ^5 of the Painleve/Gambier list [10] and the one resulting from the 
derivative of a Riccati equation. Its canonical form can be obtained by a translation 
and scaling of w. 



3. Line ARIZ ABLE MAPPINGS from the Gambier systems 



The Gambier equation is obtained as a system of two Riccati equations in cascade. 
(Cascade in this context means that the first equation contains only one variable, while 
the second one contains both. To solve the system one is thus led to integrate the two 
equations in that order, substituting the solution of the first one into the second one). 
In a discrete context the Gambier system was studied in detail (first in [6] and in full 
detail in [7]). The approach is the same as in the continuous case. One introduces two 
discrete Riccati equations in cascade: 

by + c 
ay + d 

exy + (x + r]y + 9 

x = 7, F ( 3 - 2 ) 

axy + px + 72/ + o 

with the obvious assumption that ac — db ^ 0. Eliminating y between (3.1) and (3.2) 
one obtains a single 3-point mapping for x. In this section we shall investigate the cases 
of the Gambier mapping that contain the linearizable mapping. To this end we shall 
first reduce the Gambier equation to the form (2.5). In practice this means that the 
terms quadratic in x must be absent. The conditions for the mapping to be linear in all 
x, x, x read: 

P(a( + be) = a(ce + d() 
5(a( + be) = 7(ce + d() 

a(ca + d/3) = /3(a/3 + ba) (3.3) 

7 (ecu + dp 1 ) = d(a(3 + ba) 

The solution of system (3.3) is long and tedious. We shall not enter into all the details 
but work out only the generic case and for the remaining cases give just the results. 
The generic case is based on the assumption that none of the terms appearing in (3.3) 
vanishes. By dividing the first two (or the last two) equations we obtain ad = f3^f. 
We then form the product of the second and third equation, use the previously derived 
relation and expanding we find a( = (3e. It suffices then to substitute into any of the 
four equations in order to get a final relation. The nongeneric cases are obtained by 
assuming that specific terms do vanish. Five cases can be distinguished, finally: 
a $ a55_ + bd'y — d5"f 



1) 
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e(3 
a ' 


2) 
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= P = 


3) 
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= P = 


4) 


a 




5) 


s 





0, 6 



_ 7(ce + dQ 
a( + be 



e p (3{aj3 + ba) — dj3a 
^ ~ "57' c — aa 

Once these basic conditions are obtained we still have full homographic freedom 
which will allow us to reduce further (3.1) and (3.2). Let us work out in detail the case 



corresponding to constraint 1) above. First we can translate y so as to bring 5 to zero. 
The direct consequence of5 = 0is( = /3 = c = 0. We can further translate x so as to 
get 7 = and by division put b = 1, a = 1. Finally by scaling x we can take e = 1, 
unless e = and by scaling y we can put 9=1. In the case e = we can use the scaling 
of x to put i] = 1 (unless, of course, t] = 0) and then put 9 = 1. As a result of these 
transformations case 1) can be reduced to three cases: 

y = -=-; (3-4) 

ay + d 

combined with one of the three equation for x below: 

_ xy + r]y + l 

x = (3.5) 

xy 

or 

x = (3.6) 
xy 

or 

x = — (3.7) 

xy 

Concerning the remaining cases we just give the final results. Case 2) leads to the 
following mapping: 

by - 1 

V = -^—r (3-8) 
ay + 1 

x = x y + x + (3.9) 

y 

which is the main type and also to three simpler cases 

by 

y = ~t (3.io) 

ay + l 
- xy + 1 

x = (3.11) 

y 

by + c 

y = (3.12) 
y 



and also: 



coupled to either: 

x = X -±^ (3.13) 

y 

or 

x = - (3.14) 

y 

Case 3) leads to cases identical to some of the ones identified in case 2) plus one trivial 
case corresponding to a purely linear mapping y = by + c, x = x + y. 
Case 4) yields the mapping: 



y = by + c (3.15) 

x = (3.16) 

x 

Finally case 5) is identical to case 1). We can readily remark that the nongeneric 
mapping (3.4) or (3.10) is in fact linear for 1/y. Thus the system (3.10), (3.11) can 
be readily discarded as been linear. Moreover (3.15), (3.16) is a subcase of (3.4), (3.5) 
(once y is inverted). We can thus summarize the cases to be studied in detail in the 
following list (where we have inverted y in (3.4)). We obtain just three cases, where 
a = or 1. 

I) 

y = a + dy, xx = ax + y (3-17) 

(obtained from case 1) up to a translation of y). 
II) 

y = b+-, x = - + a (3.18) 

y y 

(obtained from the simplified case 2)) and 
III) 

ty- 1 _ xy + x + ad 

y = — — 7, x = 3.19 

ay + 1 y 

(obtained from the main type of case 2)). The cases with a = are trivially integrable: 
once y is given the integration for x is obtained through a multiplicative equation of the 
form xx = F(n) or x/x = F(n), which can be linearized by just taking the logarithm 
of both sides. 

In order to investigate the integrability of the a = 1 cases we shall apply the 
singularity confinement criterion. In the case of the mapping (3.17) a singularity appears 
when y = in which case x does not depend on the value of x. No way exists for the 
mapping to retrieve this lost degree of freedom through some indeterminate form. Thus 
(3.17) for a = 1 is not integrable. We consider now the mapping (3.18) for a = 1: 

y = b+- (3.20a) 

y 

x=- + l (3.206) 

y 

The sequence of singular values of (3.20a) is y = 0, y = oo and finite values for y 
and beyond. The fact that y = induces a divergence on x, and at leading order we 
have x oc x/y. At the next step we must take into account that y oc c/y and thus the 
computation of x leads to x = 1 + x/c that is a finite value depending on the initial 
condition x. Thus the singularity is confined and the mapping (3.20) is expected to be 
integrable. Finally, we examine mapping (3.19) for a = 1. The only singularity sequence 
is y = 1/6, y = 0, y = — 1 etc. Starting with a regular x we find a finite value for x, 



while x diverges. However a close examination of the singularity shows that x is finite 
and the singularity is confined without any constraint. (Moreover, the special value 
a = 1, which leads to a more complicated singularity pattern for y, also has confined 
singularities for x). Thus the mapping (3.19) is expected to be integrable in general. 
Before proceeding further let us transform (3.19) to a form close to that of (3.20). It is 
possible to show that by translating x, x — > x — 6, and a homographic transformation 
on y, y — > y/(l — y), we can bring (3.19) to the form: 

y = A + — (3.21a) 

y 

x = - + p (3.216) 

y 

where A = (6+ l)/(a + 6), B = -l/(a + b) and p = 9-9. 

We turn now to the integration of mappings (3.20) and (3.21). Instead of restrict- 
ing ourselves to the particular forms of (3.20) and (3.21) we shall generalize our setting 
a little since this will lead to interesting results. Let us start with the more general 
Gambier system: 

y = a+- (3.22a) 

y 

x = X - + P -^ (3.226) 
y ry + s 

The application of singularity confinement follows exactly the same steps as for (3.20-21) 
and shows that (3.22) satisfies this integrabilty requirement. Thus (3.22) is expected to 
be integrable. Next we can ask for which values of p, g, r, s is the 3-point mapping for 
x is of the linearizable form (2.5). It turns out that this is possible only if the equation 
for x of the form: 

x 9 — (by , 
x=- + — 3.23 

y y 

If (h = 9 then it is possible to bring them both to zero (by tranlation of x), and the 
mapping becomes trivial. Otherwise it is possible, through a translation, to bring either 
(b or 9 to zero and scale the remaining one to unity. Thus the only interesting forms of 
(3.23) are: 

x = — (3.24) 

y 

or (with £ = x + 1) 

£=- + 1 (3.25) 

y 

One introduces the transformation x = —(1 + ujX) -1 (or equivalently £ = ujX(1 + 
u)X)~ l ) where uJulj_ = —6 and reduce the 3-point mapping resulting from the elimination 
of y between (3.22a) and (3.24) to: 

XXX + XX^^ + Xu(l-y) + l = (3.26) 

LO 



We remark that (3.26) is the generic linearizable mapping presented in its fully reduced 
form in (2.21). Thus the Gambier mapping leads again to the linearizable one obtained 
from the projective discrete Riccati system is perfect analogy to the continuous case. 

4. Conclusion 

In the previous sections we have investigated 3-point mappings that are integrable 
through linearization. Our analysis was guided by the analogy with the continuous 
situation and results of ours on N = 3 projective systems and the Gambier equation. 
We have presented a exhaustive analysis of the linearizable mapping and identified all 
its integrable forms. We thus have shown that the discrete equation: 

XXX + aXX + /3X + l = (4.1) 

can be reduced to a linear 3x3 system. The same equation, with the transformation 
X = c(l + 1/x) can be reduced to the Gambier system: 

y = a + - (4.2) 

y 

x + l 
x = 

y 

(where a, 6, c are related to the a, /?), which provides a different method for its solu- 
tion. Nongeneric cases of the general linearizable mapping were also identified and in 
particular the equation: 

x(<yx + x + 1) +7X = (4.3) 

that is the discrete derivative of a homographic mapping and constitute thus the dis- 
cretisation of equation 7^5 on the Painleve/Gambier classification. The extension of the 
present results to higher order (third and beyond) systems could be most interesting. 
We intend to return to this question in some future work. 
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